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Abstract 



^ ' We study nuclear matter and finite nuclei with a chiral Lagrangian which gener- 

H 



alizes the linear a model and also accounts for the QCD trace anomaly by means 
of terms which involve the a and tt fields as well as the glueball field (p. The 
form of the scale invariant term leading to an omega meson mass, after symmetry 
breaking, could involve coupling to (j)'^ or or some linear combination thereof. 
In fact an uj ^j^lo^^ (j)^ form is strongly favored by the bulk properties of nuclei, which 
also rather strongly constrain the other parameters. A reasonable description of 
the closed shell nuclei oxygen, calcium and lead can be achieved, although the 
spin-orbit splittings are somewhat smaller than the experimental values. 



1 Introduction 



We are still far from being able to transform the Lagrangian of quantum 
cliromodynamics (QCD) into an effective Lagrangian involving mesons and 
baryons which could be used for nuclear matter or finite nuclei. Nevertheless 
it is possible to incorporate the broken global chiral and scale symmetries 
of QCD into an effective Lagrangian. In particular the notion of (broken) 
scale invariance has been modeled by a scalar glueball potential An 
early attempt to introduce this into an appropriately modified Walecka-type 
Lagrangian was made in ref. 0. An alternative approach is to include the 
glueball potential in the Lagrangian of the linear sigma model (supplemented 
by a contribution from the vector-isoscalar omega field) [0, The latter 
is particularly attractive since it incorporates (spontaneously broken) chiral 
symmetry, another feature suggested by QCD. A common problem associated 
with these models is that the compression modulus for equilibrium nuclear 
matter is at least a factor of two larger than current estimates in the range 
200-300 MeV. Of course one can arbitarily correct this by adding terms 
to the Lagrangian, for instance and o"^, and choosing the parameters 
to produce the desired result 0. We have recently suggested a more 
satisfactory approach, however, where the form of the potential which breaks 
scale invariance is modified in a reasonable way to include a contribution from 
the cr-field as well as the glueball field, 0. This was found to yield K ^ 390 
MeV, which, while still on the large side, is decidedly more reasonable. 

There still remains uncertainty concerning the form of the omega mass 
term, im^cu^cu'^, where the explicitly dimensionful parameter must be 
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replaced by the square of a scalar field (multiplied by a dimensionless coupling 
constant) in order to satisfy the dictates of scale invariance. The question 
arises as to whether the coupling be to cr^ or 0^ or some linear combination of 
the two. In refs. Q the first choice was made. One of the purposes of the 
present paper is to examine this choice in nuclear matter and finite nuclei; we 
shall find that the properties of finite nuclei are sensitive to this question. A 
more basic purpose of this investigation is to determine whether our effective 
Lagrangian, which saturates nuclear matter in the mean field approximation, 
is able to give a sensible account of finite nuclei. This is particularly pressing 
since Furnstahl and Serot [0 have recently suggested that chiral models are 
inherently unable to give a satisfactory description of nuclei; this work was 
based on the approach of Boguta P] and did not include the glueball field, 
so there was no attempt to incorporate broken scale invariance. 

The effective Lagrangian that we employ, together with the necessary 
formalism for finite nuclei and nuclear matter is given in sec. 2. Our results 
are displayed in sec. 3, first for nuclear matter and then for the closed 
shell nuclei ^^O, ^°Ca and ^°^Pb. Our conclusions are given in sec. 4. It 
is of interest to compare the low density expansion to that of the standard 
Walecka model 0, this we discuss briefly in the Appendix. This suggests 
a simple approximation in which the glueball field is frozen at its vacuum 
value and the formalism is outlined in the Appendix. Examples are given in 
the text which show that the frozen glueball model reproduces our complete 
results quite accurately. 
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2 Theory 



We write the total effective Lagrangian in the form 

Aff ^Co-Vg, (1) 

where in this schematic equation we have separated out the scale invariant 
part, £0, from the potential Vq which induces the breaking of scale and 
chiral invariance. (The effects of explicit chiral symmetry breaking arising 
from the non-vanishing of current quark masses in the QCD Lagrangian are 
neglected.) We write Cq in terms of the chiral-invariant combination of sigma 
and pi fields, a and tt, the glueball field 0, the field of the omega vector meson 
ujfj^ and, since we are interested in finite nuclei, the vector-isovector rho field 
and the Maxwell field A^. Specifically 

+7V[7^(i9^ - g^uj^ - \gph^ ■ r - ie(l + Ts)A^) - g{a + it: ■ T-f^)]N . 

(2) 

Here the field strength tensors are defined in the usual way F^,^ = dfj,u!,^—d^u!fj,, 
^^lu — dfjhi, — 9iyb^ and /^i^ — d^A^ — d^A^. As we have remarked in the 
introduction, the form of the scale invariant uo mass term is not a priori 
obvious, so we have chosen a linear combination cx^ and (1? fields (of course 
more exotic choices are possible) and we shall examine the effect of varying 
the coefficients. In the corresponding case for the p meson we have simply 
chosen 0^, we shall comment further upon this later. 
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In order to obtain the scale breaking term Vg in the Lagrangian, we recall 
that the divergence of the scale current in QCD is given by the trace of the 
"improved" energy-momentum tensor and so the scalar potential Vg'(0, cr, tt) 
is chosen to reproduce, via Noether's theorem, the effective trace anomaly 

«;.4M*.)-i:*.||.4.,„(Ay. (3) 

where $j runs over the scalar fields {(p, cr, tt} and evac is the vacuum energy. 
The proportionality oc 0^ is suggested by the form of the QCD trace 
anomaly [pT| , 

0',{x) = , (4) 

where F^^{x) is the gluon field strength tensor and [3{g) is the usual QCD beta 
function. Eq. (3) involves a purely gluonic, color-singlet, scalar, dimension- 
four operator, and also allows for the recovery of the low energy theorems 
which follow from broken scale invariance |T^ . Then in ref . we suggested 
the form 



a, tt) = In ^ - i - 



+ li?5CV(a^ + ^^-|^) ,(5) 

where C = ~^ ^ind the subscript indicates the vacuum value. Here the 
logarithmic terms contribute to the trace anomaly and are such that eq. (3) 
is satisfied with evac = ~|-B0q(1 — 6). This requirement uniquely specifies 
the second term on the right in eq. (5). The third term is needed to ensure 
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that in the vacuum (j) = (f)Q, a = ao and tt = 0. To retain the physically 
necessary feature that evac < 0, given that B > 0, we require that S < 1. 
Provided 0<5<1,Vg'-^oo for cr or — ^> oo as is physically sensible. 
For further insight regarding the parameter 6, we recall that the QCD beta 
function with Nc colors and nf flavors is given at the one loop level by 



where the first number in parentheses arises from the (antiscreening) self- 
interaction of the gluons and the second, proportional touf, is the (screening) 
contribution of quark pairs. Recalling eq. (4), a value of 5 = 4/33 is suggested 
for the present case with Uf = 2 and = 3. Since one cannot rely on the 
one-loop estimate for (3{g), we shall examine the effect of modifying the value 



The careful reader will have noticed that Ccs, given by eqs. (2) and (5), 
does not contain the standard potential term of the linear sigma model which 
we have previously written [0, |^ in scale invariant form as fo"^ + tt^ — . 



We found that very small values of A were favored by the predicted com- 
pression modulus in nuclear matter. We have also found that departures 
from small values yield binding energies of nuclei which are much too low. 
Therefore in this presentation we make the welcome simplification of setting 
A = 0, thus eliminating this term. Note that while the term is present in 
the standard linear sigma model, it is not necessary for the breaking of chiral 
invariance. This can be carried out just as well by Vg in eq. (5). 




(6) 



of (5. 
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It is economical to write the sigma and glueball fields in terms of their 

ratios to the vacuum values, viz. 

'(r) _ a(r) 



Xir) 



(7) 



For r — s> oo we obtain the vacuum values x = ^ = ^- Now, in the vacuum, 
we require that the rho and the omega masses take their physical values so 
that we can write the mass terms in the form 



^m^x^b^ • b'^ and 



2"''uj 



(8) 



R^v^ + (1 - R^)x^ 

Thus, for the omega, taking R^^ to be 1 gives a pure z/ coupling, whereas 
the value yields a pure x coupling. The effective mass of the nucleon, in 
the usual way , is M*(r) = Mv{r). The field equations are obtained from 
Lagrange's equations and of course tt = in the mean field approximation. 
We specialise to the time independent, spherically symmetric case. The Dirac 
equation for the nucleon and the equation for the photon field are of the form 
given by Horowitz and Serot and do not need to be repeated here. As 
usual for the vector fields only the time-like components survive and also just 
the isovector z-component for the p field. The equations for the 0, a, uj and 
p can be written 

4>lDx - 2Bo{2 - 6)x + 2Bo6u = ABo[x'{\nx -6\nu)-x + 6u] 
+B,6[xiiy^ - x') + 2(X - ly)] - mlil - R^^lx - rnXx , 



a^Du - 2Bq5v + 2Bq5x = Mp^ - Bq5 



-tuIR^UqU , 



DujQ - mlujo = -g^pB + m^[i?^(z/^ - 1) + (1 - - l)]uJo 



Dbo - mlbo = -Qpps - mj(l - x^)bo , 



(9) 



where the densities ps = {NN), pB = (N'y'^N) and pa = ^{N'y^ToN) can 
be expressed in terms of the components of the nucleon Dirac spinors in the 
usual way [|^. In eq. (9) we have made the definitions -D = + f ^ and 
Bq = B(f)Q. The terms linear in the fields, i.e., the kinetic energy and mass 
terms, have been separated out on the left of these equations. The equations 
of ref. for nuclear matter are regained by setting the derivatives to zero 
so that the results are independent of 0o and (Tq, however these quantities 
are needed for finite nuclei. They are also needed to obtain the vacuum 
scalar masses and we note that the mass matrix is not diagonal. In order 
to solve the equations (9) by an iterative Green's function technique [0, 



it is necessary to go to a representation which is diagonal in the limit that 
r —>■ oo. We solve in terms of the quantities (1 — x) and (1 — z/) which go to 
zero in the same limit. 

The energy-momentum tensor can be used to obtain the total energy of 
the system in the standard way 0. Subtracting constants so that the energy 
is measured relative to the vacuum, we obtain 

Q^*^ roc f 

E = J2 ea(2jQ + 1) - 27r drr^lMups + q^ujqPb + Qphpz 

+2B,[x\\nx - 51nz/ + i) - i] + \B,5[x\2v^ - 3x') + 1] 

-mllR^u' + (1 - i?Jx>o - myblj ■ (10) 

In the first term on the right the are the Dirac single particle energies and 
ja is the total angular momentum of the single particle state. In nuclear mat- 
ter this term becomes 4 (s'w'-^o + Vk"^ + M*^^ . Making this replacement 
and using eq. (9) the expression for the energy density of nuclear matter 
given in ref. [0] can be obtained. 



3 Results 



3.1 Nuclear Matter 

In nuclear matter we have four parameters to consider. We will choose Bq 
and = ^'^2 to fit the saturation properties. The binding energy/nucleon 
we take to be 16 MeV. In refs. [0, |^ we took a saturation density of 0.16 
fm~^, but this yields central densities in lead which are too high. We have 
therefore used 0.148 fm~'^ [kp = 1-3 fm~^), which is the value favored by 



Serot and coworkers in mean field calculations |13, 0. Having fixed these 
two parameters, we then have the parameter of the uj mass term, i?^, and 
the coefficient of the quark contribution to the trace anomaly, 6, which can 
be varied. 

We show in table 1 several parameter sets that we have considered (the 
designation F here indicates the frozen glueball model where x is kept equal 
to 1; see discussion in the Appendix). One might hope that the parameters 
would show a reasonable correspondence with the independently determined 
quantities, although precise agreement would not be expected for an effec- 
tive theory such as this. In this sense the values for in table 1 are in 
acceptable agreement with the value 103 ± 36 deduced from data on the lo 
coupling constant WM and the known mass. Turning to |evac|, QCD sum rule 



studies suggest a value of (240 MeV)"^ which would not favor the extreme 



variations of 6 (sets IV and VII), although bag model estimates |16| range 
as low as (146 MeV)^. The values of Cp, given in table 1 for later use, have 
been obtained by fitting to a symmetry energy of 35 MeV. Since exchange 
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effects are known to make a significant contribution to the symmetry energy 
we shall not make a direct comparison to experiment. 

Turning to the derived quantities, nuclei generally favor effective masses 
below 0.7 [|^, which is the case for sets VI and VII. However, here the com- 
pression moduli, are in excess of 400 MeV. Pearson [|l^ has made the 
point that the leptodermous expansion, which is one method used to fit the 
data, permits values of K up to about 400 MeV with little difference in the 

fits. Defining the third derivative of the binding energy per particle as 
^ ~ ^F'i^ (^) (evaluated at equilibrium), the data indicate a linear relation 
between S/K and K (see ref. We have listed values of S/K in table 

1. Cases VI and VII are well off the allowed error band, but the remaining 
cases are quite close to, although slightly off, the band. 

The parameter sets designated with an F in table 1 correspond to the 
frozen glueball model (x = 1) and we see that they give a close correspon- 
dence with the complete results. The point is further made in figs. 1 and 
2 where the agreement is seen to be particularly good at low density, as we 
argue in the Appendix. The lower panel of fig. 2 shows that in the complete 
model X remains very close to unity over the whole region of densities con- 
sidered. We remark that the curves for u will increase linearly with kp for 
sufficiently large density (a discussion of the high density behavior was given 
in ref. 0), although it is debatable whether this theory is applicable in that 
region. 
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3.2 Finite Nuclei 

When we turn to finite nuclei two extra parameters are involved, namely 
the vacuum values of the scalar fields, 0o and cxq. Thus we have a total of 
four free parameters. Let us first eliminate two of these. We shall discuss 
changes relative to parameter set I of table 1 with cq = 110 MeV, which 
will turn out to be our preferred parameter set. First consider reducing S 
by using sets IV and V of table 1 (and reasonable values for the parameters 
not specified). We find little sensitivity in the predicted properties of nuclei. 
However, if we increase S using sets VI and VII, for which K is large, we find 
a significant reduction in the binding energy, particularly for set VII (^°Ca 
is less bound by 2j MeV, for example). We conclude that a value of 6 in 
the neighborhood of 4/33 is reasonable and fix on this value, although nuclei 
are not sensitive to the precise figure. Next consider the ratio of the vacuum 
scalar fields, ( = (Pq/ctq. Using set I we have varied ( in the range 0.7-2.1, 
which corresponds to varying the higher scalar mass, m>, between 3 and 1 
GeV. Very little change is observed in the nuclear properties. The salient 
point is that the mixing between the glueball and the sigma is small, < 3%, 
and it is the sigma which directly couples to nucleons. Henceforth we choose 
( such that m> = 1.5 GeV because this seems reasonable in view of QCD 
sum rule estimates for a scalar glueball mass in the range 1-2 GeV [1^ . The 
values of ( are listed in table 2. 

We are then left to consider the variation of R^j and ctq: a few comments 
on the expected magnitude of ctq are in order here. In the original Gell- 
Mann-Levy linear sigma model [^, a calculation of the axial vector current 
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matrix element responsible for pion decay led to the identification o"o = /tt = 
93 MeV. In a more general chiral effective Lagrangian incorporating vector 
mesons (ref. for an early review) this identification no longer follows due 



to the necessity of including the Oi axial vector meson, together with the p. 
Briefly, cross-terms of the form a^* ■ d^iz appear in the full Lagrangian which 
are eliminated by the replacement — + ^S^tt, with an appropriate 
parameter ^. This, in turn, leads to a change in the coefficient of the pion 

kinetic term, O^tz ■ d^ir, and a rescaling of the pion field is necessary to 

_i 

bring this to canonical form. One defines a renormalized field tt/j = Ztt ^tt. 

1 

Now a calculation of the full axial current yields ctq = f^^. The precise 
form of is model dependent. In the simplest model, which does not 



include the physics of broken scale invariance one finds [|2l|, using the KSRF 
relation, that = ma^/nip = \p2. We do not expect either this precise 
functional form or value to be maintained in more realistic models, but we do 
anticipate Zn > I. Accordingly, do is kept as a free parameter, which we can 
presume to be larger than 93 MeV. Two further remarks are in order. First, 
the additional terms involving the ai will not contribute in the mean-field 
approximation. Second, from eq. (2) after symmetry breaking, the nucleon 
mass is given by M = gao- The physical pion- nucleon coupling is defined in 
terms of the renormalized pion field to be g^NN = Z^g. So, in fact, one still 
has M = gao = g-^NNf-K as the approximate form of the Goldberger-Treiman 
relation. 

With this in mind, we first consider variation of ctq keeping R^^ = 0. The 
results are given in table 2, where the notation 1/93, for example, implies use 
of parameter set I from table 1 with ctq = 93 MeV. The theoretical values of 
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the binding energy/particle include a correction for the cm. kinetic energy 



[^ . The charge densities, and therefore the radii quoted, are corrected for 
the finite size of the proton [|13[ and for cm. effects [^]. We see that the 
choice (Jo = 93 MeV (fourth row in table 2) yields poor binding energies, radii 
which are too small and, as shown in fig. 3 for ^°Ca, an oscillatory charge 
density. This figure indicates that increasing ctq to 110 MeV improves matters 
greatly. One can view this as reducing the magnitude of the derivative of v 
in eq. (9) or as reducing m< to ~ 500 MeV (Furnstahl and Serot 0] stress 
that values in this neighborhood are to be preferred). In fact as regards 
the charge radii and densities the value 120 MeV gives even better results, 
however we see from table 2 that in this case the binding energies/particle 
are l-l| MeV too low. Therefore we prefer the parameter set I/llO. 

Now we consider variation of R^^ with ctq = 110 MeV. As R^^ is increased 
from zero to 0.5 or 1.0 we see that the binding energies become too large and 
shows the wrong A-dependence and also the radii become too small. Table 2 
also shows that the mass of the lighter scalar meson, m<, increases sharply. 
Further, as we see for ^osp]^ ^ charge density begins to develop 

oscillations which are not present in the data. Similar problems were noted 
by Furnstahl and Serot ||^ who used R^ values of 0.76 or 1, although their 
model differs significantly from ours. As pointed out by Price, Shepard and 



McNeil this is likely to signal the onset of a situation where the ground 
state of nuclear matter is no longer uniform, but exhibits periodic density 
fluctuations. We conclude that nuclei strongly favor R^ = 0. 

We have also compared the case I/llO with the frozen glueball approxi- 
mation, case IF/110, in table 2. There is a close correspondence between the 

13 



results, as there is for the predicted charge densities which we illustrate for 
1^0 in fig. 5. 

Finally we turn to the single particle energies. The levels near the Fermi 
surface in lead are shown in figs. 6 and 7 for neutrons and protons respec- 
tively. The relative positions of the neutron and proton levels are sensitive to 
the treatment of the p meson JTSf . We have used a scale- invariant 0^ coupling 
in the mass term, but the results differ negligibly if a standard mass term 
is used. However there is a difference if a cr^ coupling is used because this 
reduces and leads to ~ 2 MeV less (more) binding for protons (neutrons). 
Since figs. 6 and 7 clearly show that this is undesirable we have not consid- 
ered this possibility further. The results with the sets 11/110 {Ruj = 0.5) and 
1/93 show immediate problems since the major shell closure is incorrectly 
given due to the position of the proton 1/19/2 and the neutron liii/2 levels. 
This problem disappears with the preferred parameter set, I/llO, and the 
single particle spectrum is much more reasonable, although a larger proton 
1/^9/2 — 3s 1/2 splitting would be desirable. This may reflect the spin-orbit 
splittings which are, on average, 65% of the experimental values here and 
also in the other nuclei we have considered. They could be improved by 
reducing the effective masses in table 1 below 0.7. This can be achieved by 
taking a small negative value for i?^, although we find that this leads to a 
larger compression modulus and lower binding energies. 
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4 Conclusions 



Our primary purpose has been to see whether an effective Lagrangian which 
includes the brealcing of scale and chiral invariance can make sensible predic- 
tions for finite nuclei. Finite nuclei are expected, and found, to be a much 
more severe test than nuclear matter. In this endeavour the form of the scale 
invariant mass term for the uo meson was found to be critical and a pure cou- 
pling to the glueball field, cfp'uj^uj^, was required, although a small effect from 
the sigma field cannot be ruled out. A similar coupling for the p mass term 
was also favored. It is interesting that the properties of nuclei can be used 
to provide such strong constraints on the form of the couplings of the scalar 
glueball to other mesonic fields appearing in the chiral effective Lagrangian. 

It was also found to be necessary to increase the value of the vacuum value 
of the sigma field, ctq, from the naive expectation of 93 MeV to somewhere in 
the neighborhood of 110 MeV. With these caveats the bulk properties of O, 
Ca and Pb were quite well accounted for and, while the spin-orbit splittings 
are about 65% of the desired values, the basic structure of the single particle 
spectra was reasonable. As a by-product we were able to define a simplified, 
frozen glueball model which was able to reproduce the results of the complete 
model quite accurately, but, by itself, would seem highly arbitary. 

It is natural to compare this approach to the standard Walecka model ex- 
tended to include non-linear and terms (as for example in refs. ^). 
The agreement is better than found here, particularly as regards the single 
particle energies. Nevertheless, it is remarkable that the present effective 
chiral Lagrangian, whose form is motivated by QCD, can give an acceptable 
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description of both nuclear matter and finite nuclei in the mean field ap- 
proximation. It remains to be seen whether fine tuning can further improve 
the agreement here. For example, a term of the form (cu^cj'^)^ is the sim- 
plest among many possibilities, consistent with scale invariance, which could 
be added to the Lagrangian. We are currently investigating the effects of 
including such a term. 

After this work was completed we received a preprint from Furnstahl and 



Serot l2§] who have studied the same problem, restricting themselves to an 
uo mass term of the form a'^u^u'^. In agreement with our conclusions, they 
point out that the results for nuclei are poor. As we have stressed, nuclei 
require that the coupling be to the glueball field rather than the sigma field. 

We thank G.F. Bertsch and A. Vischer for useful comments and C.J. 
Horowitz for giving us a copy of his computer code for finite nuclei which 
was suitably modified for present purposes. We acknowledge partial support 
from the Department of Energy under contracts No. DE-FG02-87ER40328 
and DE-AC02-83ER40105. A grant for computing time from the Minnesota 
Supercomputer Institute is gratefully acknowledged. 
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Appendix 

Low Density Expansion and Frozen Glueball Model 



It is of interest to study tlie low density expansion of eqs. (9) and (10) 
in the case of nuclear matter. Examination of eq. (9) indicates that we can 
expand 



i/ = 1 + apB{l + hkl) + 0{kl) , 

x = i + cpB(i + d4) + c(4) , 



(11) 



where pb is the baryon density and kp the Fermi momentum. Further to 
this order the coupling between cuq and v and x can be neglected. Then the 
equation of motion for % indicates that 

c 5 
a 



2-(5 ' 

The equation of motion for i/ yields 

M(2 -S) , 3 



(12) 



a — 



ABo5(l - 5) ' 10M2 ■ 

Substituting into eq. (10), we find the energy/particle 



(13) 



E 
~A 



M+ '^"^ 



3kp 



lOM 56M3 
M\2 - S) 



+ 



qIpb 



2m2 



Pb\1 



3k'p 



+ 



kp' 
M 



(14) 



85o5(l-5)^"V 5M2^ 
The first three terms are the energy of a Fermi gas and the succeeding terms 
are the cu and scalar meson contributions. This is of exactly the same form 
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as the standard Walecka model |p, except that the scalar meson mass is 
replaced by 

^ = — 277; =m^my ~ m< . (15) 

V ^ol^-^J / [2Bo{2-6)]2 

where m< {m^) are the lower (upper) eigenvalues of the vacuum mass ma- 
trix (see for an explicit discussion). The mass, m, is approximately m<, 
provided that (po/cxo is not too large so that the mixing of the sigma and 
glueball states is small. 

Now the fact that d = b in eq. (12), means that, to the order we are 
considering, we can obtain the correct result for u with x = 1) provided that 
Bq is replaced by 

= (16) 

In this approximation, we see that the scalar potential, Vg, takes the simple 
form 

Vg(0, a, TT = 0) - e^ac = |5o5(z/2 - 1 - In z/2) . (17) 

In eq. (9), with x = 1; the glueball equation of motion can be dropped and 
we are left with the single scalar equation 

alDu - 2Bo6iy = Mps - Bq5 ^^ ^ - mlR^uolv . (18) 

Taking care that the energy is consistently obtained we find 

E = Y^ ea{2ja + 1) - 27r / drr'^l Mvp^ + g^jUJopB + fi'p&oPa 
a Joy 

+2Bo5\n V - mlR^v'u?^ . (19) 
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Using the low density expansion for nuclear matter and replacing Bq by B'q 
in eqs. (18) and (19) this expression gives the same result as the complete 
result of eq. (10). In practice to saturate nuclear matter correctly a small 
adjustment of is needed and Bq is only approximately given by (16). 
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Table 1 

Values of the parameters and the derived quantities for nuclear matter 



Set 






1 

1 Cvac 1 

(MeV J 






M 


(iVleV J 


s 

K 


1 


U.U 


A 

4 


ZOO 


lo6 


116 


0. / 1 


QOQ 

ooo 


6.0 


IF 


0.0 


4 


231 


152 


123 


0.72 


341 


5.6 


II 


0.5 


4 


255 


82.6 


127 


0.80 


356 


5.9 


III 


1.0 


4 


269 


51.3 


132 


0.84 


377 


6.1 


IIIF 


1.0 


4 


265 


49.2 


135 


0.85 


356 


5.7 


IV 


0.0 


1 


336 


153 


122 


0.72 


350 


5.8 


V 


0.0 


2 


281 


154 


119 


0.72 


360 


6.0 


VI 


0.0 


8 


194 


158 


107 


0.69 


442 


7.5 


VII 


0.0 


16 


160 


156 


88 


0.65 


666 


10.6 
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Table 2 

Bulk properties of nuclei for various parameter sets 









r\ 
U 




L^a 




TDK 

r D 




Set 


m< 


c 


BE 
A 


fch 


BE 
A 


fch 


BE 
A 


fch 




(MeV) 




(MeV) 


fm 


(MeV) 


fm 


(MeV) 


fm 


Experiment 


7.98 


2.73 


8.45 


3.48 


7.86 


5.50 


I/llO 


506 


1.4 


7.35 


2.64 


7.96 


3.41 


7.33 


5.49 


lF/110 


508 




7.86 


2.62 


8.35 


3.40 


7.54 


5.49 


1/93 


598 


1.4 


9.41 


2.52 


9.38 


3.32 


8.04 


5.45 


1/120 


464 


1.4 


6.25 


2.72 


7.16 


3.47 


6.92 


5.52 


11/110 


598 


1.6 


10.08 


2.49 


9.85 


3.29 


8.22 


5.45 


Ill/no 


664 


1.8 


10.98 


2.44 


10.45 


3.26 


8.48 


5.44 
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Figure Captions 



Fig. 1. The binding energy/particle as a function of density for nuclear 
matter; density p is given as the ratio to the equilibrium value po = 0.148 
fm~^. The parameter sets used are indicated and the designation F indicates 
the frozen glueball model. 

Fig. 2. The fields v and x corresponding to fig. 1. 

Fig. 3. Comparison of the experimental charge density [^] for ^°Ca with 
theoretical predictions for the parameter sets indicated. 
Fig. 4 As for fig. 3, but for ^ospb. 

Fig. 5 As for fig. 3, but for ^^O . The data are from ref. p6| . 
Fig. 6. Occupied and unoccupied neutron levels near the Fermi energy in 
^°^Pb. The experimental data are compared with predictions for the param- 
eter sets indicated. 

Fig. 7. As for fig. 6, but for protons. 
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